The Ngai model of relaxation: An alternative 
for the study of an interacting system of small 

particles 

Running Title: Relaxation in an interacting fine 
particle assembly 
Subject Classification: 75.60 
J. Castro 
Departamento de Fisica Aplicada 
Universidad de Santiago de Compostela 
E-15706, Spain 

February 1, 2008 

Abstract 

It is shown that the Ngai model of relaxation accounts for the 
main experimental results of the magnetic relaxation in an assembly 
of interacting fine particles. 
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1 Introduction 



The Ngai model for relaxation in complex systems was first introduced in 
1979 |], U- Since then, it has been found to offer an accurate description of 
a wide variety of systems ||. However, its use in magnetism is more recent. 
See reference [|| and other references therein 

The fundamental difference between the Ngai model and previous models 
of relaxation is that these earlier models treat the effect of complexity on the 
relaxation as being essentially static, in the sense that the effect of complexity 
is to change the relaxation time or to produce the presence of a distribution 
of relaxation times. In contrast to this, in the Ngai model, the relaxation in 
complex systems, is dynamical in nature. 

From this view-point, according to Ngai [f|, a relaxing complex system 
consists of three parts: 

i) Individual primary species (PS) of which are of interest for the relax- 
ation. 

ii) A heat-bath whose interaction with the PS provides a primary mech- 
anism of relaxation. 

hi) Other relaxing species, X, whose interaction with PS, the PS-X cou- 
pling, slows down the relaxation process. This is the main manifestation of 
complexity. 

In this paper we propose another application of the Ngai model to a 
magnetic system: namely that of an interacting fine-particle assembly. In 
this system, the PS are the individual particles whose coupling with the other 
particles in the assembly slows down the relaxation process. In preceding 
attempts to tackle this problem 0, U, the only effect of the presence of 
interparticle interactions is a change in the relaxation times. A discussion 
of these results in the light of the present model, together an evaluation of 
experimental results using this model, are also presented. 

In this paper we show that many of the experimental results for an inter- 
acting system of small particles come out naturally within the Ngai model. 
In particular the model can explain, in a natural way, not only the abserved 
differences for a system of interacting and non-interacting particles via the 
variation of the inverse of the blocking temperature with the logarithm of 
the measuring time, but also the values encountered for the relaxation times 
(table 1). The organization of the paper is as follows: in section 2 we review 
the Ngai model for the relaxation in a complex system. In section 3 we apply 
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this model to a system of interacting particles and show that the behaviour 
discussed previously can be accounted for by this model. Section 4 is de- 
voted to the calculus of the magnetic viscosity in two different extentions of 
the model, the case of a distribution of particle-sizes and also the case of a 
distribution of switching fields. In the former, we encounter the ubiquitous 
logarithmic relaxation law, whereas in the latter our results agree with those 
of Charap[10| and therefore the model is compatible with the experimental 



results of Oseroff et al 11 



2 The Ngai model 

After the justification of the necesity of the study of non simple-exponential 
relaxations we sumarize briefly the more important results of the Ngai model . 
For more information about the theoretical derivations and the experimental 
succeses of the model we refer to the bibliography ||. 

It was shown by Ngai et al [H] that all the permissible relaxation functions 
have to verify the Paley- Wiener theorem of the Fourier transform theory. It 
is easy to show that the simple exponential relaxation function violates this 
theorem. Sice a single, linear exponential form is unphysical, as demon- 
strated, the idea of a superposition of exponentially decaying functions may 
be also ruled out as a description of relaxation phenomena. There is therefore 
the neccessity to move on non simple-exponential relaxations. 

The Ngai model || proposes the existence of a temperature-insensitive 
crossover time, t c , separating two regimes in which the dynamics of the 
relaxation are different. The relaxation function is a linear exponential 
exp(—^) for t < t c and cross-over at t ~ t c to a stretched exponential 
form exp(— for t > t c . The coupling parameter n lies in the range 
< n < 1 . Notice the appearance of two different characteristic time-scales, 
T and r*, for the two regimes. These two time-scales are not independent, 
but continuity of the relaxation function at t = t c implies that: 

r* = (i) 

Notice an important consequence of equation [l] which will be important in the 
following, if tq describes the physics of some motion over an energy barrier 
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with activation energy E A - 



r = r QO exp(^) (2) 



Then for equation [I], r* describes a macroscopic relaxation time given by: 

r* = r* ezp(||) (3) 



where: 
and: 

Where E%> E A . 



kT' 

E!=^~ (4) 
1 — n 

r* = r^tr 1 (5) 



3 Application of the Ngai model to the re- 
laxation of a system of small interacting 
particles 

There are several excellents treatment of the magnetic properties of an en- 



semble of particles. We shall, therefore, not enter into detail here (|TJ, . 

We consider, firstly, an ensemble of particles of equal volume V. For 
simplicity we assume that the particle anisotropy is uniaxial and that the 
particle's easy axis of magnetization is parallel to the field direction. Then 



the energy of one of the particles in a field H will be [|TJ 



E = KVsin 2 9 - M s VHcos9 (6) 

where K is the anisotropy constant, M s is the bulk saturation magnetization 
and 9 measures the orientation of the magnetic moment of the particle with 
respect to the field. This energy gives rise to two energy minima separated 
by an energy barrier: 

E A = ^M S H K (1 - J-) 2 = e A V (7) 
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At a finite temperature, we have a finite probability of a thermally activated 
motion over the energy barrier. The characteristic time-scale of the motion 
will be given by the expression: 

r = r 00 exp(E A /fcT) (8) 



where 10 9 — 10 12 s [17|. For such a system, the Ngai model pre- 



dicts the existence of a cross-over time, t c , independent of temperature, to 
be determined from the experiment, separating two different regimes of the 
dynamics of the relaxation. The relaxation function is a linear exponential 
exp(—t/T ) for t < t c and has a stretched exponential form exp(— (t/r*) 1- ™) 
for t > t c , where < n < 1. The relation between r* and tq is given by 
equation 1. 

An important consequence of the model can be now pointed out if, as a 
definition of the blocking temperature, we use the usual criterion that the 
characteristic time-scales of the relaxation will be of the order of the measur- 
ing time. We then expect that the logarithm of the measuring time plotted 
as a function of 1/Tg will be a "piece-wise" straight line with two different 
slopes, !^ for t < t c and kB ^ 1 A _ n ^ for t > t c . The slope corresponding to the 
small measuring times will be smaller than that corresponds to the longer 
measuring times. This fact was encountered experimentally in a detailed 
analysis of the dynamics of a system of small iron particles, magnetically 
interacting and dispersed in an amorphous alumina matrix by Dormann et 
al ||. We reproduce here data from a figure from that paper, figure 2, which 
shows the predicted behaviour (figure 1). From this figure, we estimate that 
t c is in the range 10~ 10 — 10~ 5 s. However these results are not immediately 
comparable with the above developed theory because they correspond to 
samples with different sizes and interactions. 

A better comparison with the theory can be made with the data of J. L. 
Dormann et al||15| in a system consisting of 7 — Fe20s particles dispersed in 
polyvinilic alcohol (figure 2) Q In this case we have a series of samples with 
the same size and different interactions including the case of zero interaction 
(sample IF). The interaction increases from right to left. From the sample 
with zero interaction we can estimate « 10~ 10 s. By comparing the slopes 
of sample IF with the slopes of the other samples, we can calculate the value 



1 I am indebted to Dr. D. Fiorani for pointing this possibility out to me and for making 
available his experimental data 
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Sample 


n 


T*oo(th)(s) 




CH 


0.4 


4 x 10~ 12 


io- 12 


IN 


0.73 


3.87 x 10- 16 


10 -ie 


FLOC 


0.76 


3 x 10" 17 


io- iV 



Table 1: Comparison between the theoretical and experimentally obtained 
values for r* . See text. 



of n. Using these values of n, the afore-mentioned value of r M and choosing 
t c « 10~ 9 , we can calculate theoretically the value of r^. These values are 
compared with the experimentally obtained values in table 1. We observe 
excellent agreement between theory and experiment. 

Here we wish to stress the importance of this fact. The existence of the 
cross-over time, which separates two different time dependences with times- 
scales related by equation 5, is a stringent condition for the verification of the 
Ngai model. This is the first time that this cross-over has been encountered 
in a magnetic system. 

A further important consequence is that the Ngai theory can explain the 
discrepancy between the size of the particles measured using two different 
time scales. This will be investigated in a further publication ]16 . 



After showing that the Ngai can describe the relaxation of a system of 
small interacting particles we dedicate the remainder of this paper to the 
extention of the model to a calculation of the magnetic viscosity in the case 
of a distribution of particle sizes and a distribution of switching fields. 



4 Calculation of the magnetic viscosity 

A typical viscosity measurement in a small particle assembly runs as fol- 
lows: After saturation in a strong field, a small measuring field is applied 
in the reverse direction. The magnetization is then measured as a function 
of time. Typical times between application of the saturation and the first 
measurement are of the order of seconds and the time-scale of the measure- 
ment is of the order of minutes. In the previous section, we estimate that 
t c ~ 10~ 10 — 10~ 5 s (figure 1) and thus we expect that the relaxation function 
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Figure 1: Variation of the blocking temperature, Tg, with the logarithm of 
the measuring time for different samples of iron particles, according to [7]. 
The average particle size and the interparticle interaction strength increases 
on going from right to left. 
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Figure 2: Variation of the blocking temperature, Tg, with the logarithm of 
the measuring time for different samples of 7 — Fe20^ particles, according to 
[14]. The interparticle interaction strength increases from right to left. The 
average particle size is held constant. 
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for the present experiment will be of the stretched exponential type. Using 
the boundary conditions, M(t = 0) = — M and M(t = oo) = M the time 
dependence of the magnetization can be written in the form: 

M(t) = M - 2M exp{-(—) 1 -"} (9) 

We now extend the model to two cases: distribution of particle sizes and 
distribution of switching fields. 



4.1 Distribution in the particle size 

We now consider the fact that we have a distribution in the particle size, 
f(V) such that: 

rco 

f(V)dV = 1 (10) 



In this caes, our relaxation law transforms to: 
M(t) 



M 



roo f 

/ (l-2exp{-(-) 1 - n })f(V)dV (11) 
Jo r* 



As a simple example we consider the distribution [[18 



f(V) = I ^ Vl < V < V2 (12) 
( J \ otherwise 

For this distribution function, the relaxation function can be evaluated 
in terms of the exponential integral: 



poo p y 

E l (x) = / dy (13) 

Jx y 

A simple calculation gives: 



kT 

M(t) = M - 2M — —AE t (z 2 ) - E t ( Zl )) (14) 

ZA\V2 — Vi) 



where: 

' ' / 1.2 (15) 



\1— n 

e A Vi\) 
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We stress here that this opens up a new possibility for the analysis of the 
relaxation data. 

For Vi, V 2 , n and t such that z\ » 1 and z 2 << 1 we have: 

M(t) = M - - n )logr 2 - 0.577 - (1 - n)/o^) (16) 

eA[V 2 - Vi) 

where: 

r 2 = r^e fc s^i-«) (17) 

We have therefore encountered in this limit the ubiquitous logarithmic relax- 
ation. In the literature, the parameter S is usually defined by: 

S = dM W M ° ( 18 ) 

dlogt v ; 

which in this case is given by: 

2k B T(l-n) 

e A (V 2 - V 1 ) [U) 

Two important comments can now be made. Firstly, the value of 5* depends 
on n and therefore on the particle interactions and secondly, there are several 
possibilities to obtain the logarithmic relaxation: (i) from equation (9) and, 
assuming that (1 — n)log(t/T±) is very small, we have|l9|]: 

M(t) = M - 2M e + 2(1 - n)M elog(t/n) (20) 

and (ii) assuming a distribution of volume sizes. 

4.2 Distribution of switching fields. 



Following Charap[T(| we consider an ensemble of interacting magnetic parti- 
cles with uniaxial anisotropy, thus with a switching barrier given by equation 
7 and a switching field distribution p(Hx)- We can transform the p(Hx) in 
a density of energy barriers G(E\) through: 

f) P* 

G{E\) = {q^)- 1 p{Hk) (21) 
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In the calculation of G(E* A ) we use the Charap method. For a given Hk, 
only for fields H ~ Hk we encounter appreciable switching. In this range of 
fields we can write approximately: 

* » W-w^ < 22 > 

which simplifies the evaluation of equation 21. To take a specific case, we 
assume a switching field distribution of Lorentzian form: 

= ,A ( i + (W^> (23) 

where A measures the half-width of the distribution at the half-maximum. Then 
G(E A ) can be written as: 

GiE:,) = — ^ = (24) 

The magnetization is now described by: 
M(t) 



M 



Defining: 



we have: 



1 - 2 f° G^e - ^ 1 ""^ (25) 
Jo 



v = ^ (26) 



M(i) if(l-n) f 00 exp(-'* 



vrA io rfX ^(i + (£il±£jy^h^) 2) (28) 



M, 
and: 

gl-j; lM r ex P (-^(l-n) 2 -(4) 1 -.e^d^ ) 
lJ vrA V io ^ (1 + ( g (i + (i-n)Vx)-^ )2 

(29) 
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We have calculated numerically the dependence of S(t) for t = 100s with 
the applied field at different temperatures for several values of n. A typical 
example, for n = 0.5, is shown in figure 3. The choice of parameters is the 
same as that of Charap: r = 10" 9 s, V = 10" 15 cm 3 , H c = 200Oe, A = lOOe 
and t c = 10~ 9 s. Several comments can now be made. S max , the maximum 
value of S, decreases with n and the curve is broader that the corresponding 
curve for n = 0. We have fitted the temperature dependence of S max to a T a 
law (figure 4) and have found a = 0.443. This does not deviate significantly 
from the n = estimate a = 0.478. From a similar calculation given by 
Charap, we expect in the first approximation, a dependence of S max with T 
of Ta. Finally, we mention that these results are in qualitative accord with 
the experimental results of Oseroff el al [II] . 



5 Conclusions 

After showing that the Ngai model of relaxation can account for the main 
experimental results of the magnetic relaxation in an interacting fine particle 
assembly, we have extended the model to the case of a distribution of particles 
sizes and switching fields. We mention finally that we look forward to this 
model being used for an analysis of the magnetic relaxation in an ensemble 
of small magnetic particles. 
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Figure 3: Expected dependence of S with H at different temperatures ac- 
cording to equation 29 
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Figure 4: S max as function of temperature for n=0.5 and fitted to a T a law. 
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